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We have used a Hartree-type electron-helium potential together with a den- 
sity functional description of liquid 4 He and s He to study the explosion of 
electron bubbles submitted to a negative pressure. The critical pressure at 
which bubbles explode has been determined as a function of temperature. It 
has been found that this critical pressure is very close to the pressure at 
which liquid helium becomes globally unstable in the presence of electrons. 
It is shown that at high temperatures the capillary model overestimates the 
critical pressures. We have checked that a commonly used and rather simple 
electron-helium interaction yields results very similar to those obtained us- 
ing the more accurate Hartree-type interaction. We have estimated that the 
crossover temperature for thermal to quantum nucleation of electron bubbles 
is very low, of the order of 6 mK for 4 He. 
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1. INTRODUCTION 

It has been recognized that liquid helium is especially well suited for 
homogeneous cavitation studies. On the one hand, it can be prepared in a 
high-purity state, avoiding heterogeneous cavitation driven by impurities in 
the liquid. Besides, experimental techniques have been developed I 1 ! 2 ! 3 ! 4 ! that 
focus a short burst of ultrasound into a small volume of bulk liquid, thus 
preventing cavitation on defects at the walls of the experimental cell. On 
the other hand, helium remains liquid down to zero temperature (T). This 
allows to address quantum cavitation, a phenomenon that may appear at 
very low temperaturesP^ 

Heterogeneous cavitation produced by impurities purposely introduced 
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in the liquid is also interesting by itself,"-" and in a series of recent experiments 
the case of heterogeneous cavitation cau sed b y electrons (electron bubble 
explosions) has been addressed in detailP^^I Another interesting case of 
heterogeneous cavitation in liquid 4 He thoroughly studied is that ca used by 
the presence of quantized vortices acting as cavitation seedsP^ 12 * 13 * 

In this work we attempt a theoretical description of electron bubble 
explosions using a T-dependent density functional approach we have em- 
ployed in a series of studies on cavitation and nucleation in liquid helium 
(see Ref. El f° r a comprehensive review) in conjunction with a realistic 
electron-helium (e-He) effective potential. Our results are in agreement with 
those obtained in Refs. and for He, and in Ref. EI|for 3 He, covering a 
wider temperature range. 

This paper is organized as follows. Sections 2 and 3 are mostly de- 
voted to 4 He and to the general formalism. In Sec. 2 we present the results 
obtained using the capillary approximation. In the case of electron bub- 
bles, this approximation has been considered realistic enough to yield semi- 
quantitative results for critical pressures and thermal to quantum crossover 
temperatures, and constitutes a useful guide to the results obtained within 
density functional (DF) theory. In Sec. 3 we present the DF plus Hartree 
electron-effective potential approach together with the results obtained for 
4 He using this method. In Sec. 4 we present the results obtained for 3 He, 
and a brief summary is presented in Sec. 5. A preliminary version of part 
of this work has been presented elsewhereP^ 

2. CAPILLARY MODEL 

In this simple model, the electron is confined in an impenetrable spher- 
ical well potential of radius R. The total energy of the e-He system can be 
written function of the radius 

, s 7r 2 ?i 2 9 4 o f-le 2 

where P is the pressure applied to the system, a is the surface tension of 
the liquid, and e is its dielectric constant. The first term is the energy of 
the electron in the ground-state of an infinite well potential of radius R. 
For 4 He the last term can be evaluated taking e = 1.0588 (Ref. IT%|) and 
£ = 1.345 (Ref. 117)1 . Its effect is small and it will not be considered in 
the following. We have also checked that the effect of including a curvature 
energy term in Eq. (^TJ) is small. On the contrary, the effect of the surface 
tension on any cavitation process is crucial, and quantitative results can only 
be obtained with the use of the correct value of a. In the following, we will 
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takP a = 0.272 K A" 2 , instead of the value a = 0.257 K A 2 given in 
Ref. [2U] that we and other authors have used in the pastP^^ED \y e want 
to mention that the value of Roche et alP^ agrees well with that of Guo et 
alP3 obtained long time ago. 

When P > 0, Eq. Q has an absolute minimum located at R m in = 18.9 
A at P = 0. This configuration corresponds to a stable electron bubble. 
When the liquid is depressed below its saturation vapor pressure, the abso- 
lute minimum becomes local, and U(R) also displays a local maximum. The 
metastability region in the P — T plane extends between the liquid-vapor 
coexistence line down to the instability line. Consequently, metastable bub- 
bles can be formed at positive and negative pressures as well. In that re- 
gion, the electron bubbles are metastable, and an energy barrier of height 
AU = U (Rmax) — U(R m i n ) appears which can be overcome either by thermal 
activation above the barrier or by quantum tunnelling through it. 




P(bar) 

Fig. 1. Energy barrier height AU (K) (dashed line) as a function of P (bar) 
of an electron bubble in liquid 4 He at T = in the capillary model. Also 
shown is the corresponding DF result (solid line). The crosses along the DF 
curve have been obtained using a simpler electron-He interaction, see text. 

The height of the energy barrier AU is displayed in Fig. \T\as a function 
of P at T = K. It can be seen that if the negative pressure is large enough, 
the barrier eventually disappears and the system becomes globally unstable. 
This happens at an instability pressure P u given by the expression: 
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For the parameters we use, P u = —2.12 bar, and the radius of the correspond- 
ing electron bubble is R u = 28.2 A. Had we taken into account the last term 
in Eq. (^Q), we would have obtained P u = —2.24 bar and Ru = 27.1 A. The 
instability pressure can be compared to the spinodal pressure at which pure 
4 He liquid becomes macroscopically unstable, P sp = —9.20 bar This value 
of P sp is consistent with the value obtained by other authors! 9 * 23 * 24 * 

A necessary condition for the validity of the capillary model is that the 
metastable electron bubble has a fairly large radius. Only in this case one 
may split the system into a volume and a surface region which justifies the use 
of Eq. (^Q). This happens when the critical configurations are large enough, 
as for example in the case of cavitation occuring near the liquid-vapor co- 
existence line! 23 * 25 ] However, in the case of electron bubble explosions the 
situation is more complex because the e-He interaction is strongly repulsive 
and small changes in the electron wave function may cause a sizeable effect 
on the metastable bubble. It is obvious that the bubble configurations used 
in the capillary model have little flexibility, and consequently, the validity 
of this approximation should eventually rely on the comparison with more 
realistic methods, as the DF approach. This will be done in Sec. 3. 

In most cases, cavitating liquids undergo phase separation before they 
reach the stage of global unstability. This proceeds through the formation of 
critical bubbles either by thermal or by quantum activation. We will show 
that for electron bubbles, the critical pressure P cr at which it happens is 
very close to P u , a result already obtained in Ref. El As a general rule, 
the presence of impurities in the liquid results in a sizeable decrease of \P cr \. 
Quantized vortices in liquid 4 He play the same role as impurities, and their 
presence also decrease |P cr | (Refs. I11I12|) as well as the degree of critical 
supersaturation in 3 He- 4 He liquid mixtures.^ 

Within the capillary approximation, the dynamical evolution of the elec- 
tron bubble can be parametrized by one single collective variable, namely 
the radius of the bubble. In this case, it is rather simple to describe quantum 
and thermal cavitation regimes on the same footing, continuously passing 
from one to the other. This is accomplished by the use of the functional 
integral method (FIMj^HEZI thoroughly described in Ref. 1131 We now recall 
its essentials. 

The nucleation rate J for a thermally activated process, i.e., the number 
of critical bubbles formed in the system per unit time and volume, is 

J T = JoTexp(-AU/k B T) , (3) 

where ks is the Boltzmann constant. The prefactor Jqt depends on the 
dynamics of the process, and it is of the order of the number of cavitation 
sites per unit volume (the number of electrons per unit volume, n e , in the 
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present case) times an attempt frequency vt- For simplicity, we take vt = 
ksT/h, where h is the Planck constant. 

At low enough T, thermal activation is no longer possible. However, 
cavitation may proceed by quantum tunnelling. The transition from one 
regime to the other is very abrupt, so that a thermal-to-quantum crossover 
temperature T* may be defined by indicating whether nucleation takes place 
thermally (T > T*) or quantically (T < T*). In the limit of zero temper- 
ature, cavitation is purely quantal, but for T* > T > thermally assisted 
quantum cavitation is the physical process. 

For T < T* the tunnelling rate is 

J Q = J 0Q eMS Q /h) , (4) 

where V = exp(—S^/h) is the tunnelling probability, is the quantum 
action, and the prefactor Jqq is again of the order of the number of nucleation 
sites per unit volume times an attempt frequency vq which can be estimated 
from the zero point motion of the system about the metastable equilibrium 

position Rmin- 

An analytical expression for T* is obtained which involves the second 
derivative of U(R) and the value of the collective mass of the bubble M(R) 
at the maximum of the cavitation barrier, R ma x 



k B T* 



1 -1 d 2 U 




M(R max ) dR? 


Rmax 



(5) 



If the motion is irrotational and the fluid incompressible, the collective mass 
M depends on R afPIED 

47T 

M(R) = -^R 3 m HePb , (6) 
n 

where pi, is the particle density of the metastable bulk liquid and m# e is the 
mass of a helium atom. Eq. (J5J) shows that the value of T* is determined by 
small variations around R ma x- An expression for the attempt frequency in 
the quantum regime can be also worked out easily: 



- J_ / 1 d 2 U 
VQ ~ 2^V M{R min ) dR 2 



(7) 

■R-rnin 



Using Eq. (JSJ) we have determined T*{P) and the results are shown in 
Fig. [21 This figure shows that in the capillary approximation, irrespective 
of the pressure, above T ~ 10 mK the cavitation process is thermal and 
not quantal. To determine which T* corresponds to the actual experimental 
conditions, one has to look for the intersection of the curve T*(P) with the 
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line that results imposing that critical bubbles nucleate with appreciable 
probability inside the experimental volume during the experimental time: 



1 



texpVexp 



—S 1% 



-AU(P)/k B T* 



(8) 



where t exp and V exp are the experimental time and volume, respectively, and 



Jqq = n e UQ. Taking t exp ~ 10 



Vexp ~ 10 



5 cm 3 , 



10 6 



cm 



and 



the value vq ~ 10 9 s _1 obtained from Eq. ( |7| ). one gets T* = 4.7 mK at a 
pressure slightly above P u . At T* = 4.7 mK one has S® = 11.5 K, so that 
the use of the formalism of Ref. well justified Had we used the WKB 
approximation^ to estimate T* , we would have obtained similar results. In 
particular, the maximum of T^ KB is 9.6 mK. We will see in Sec. 3 how 
these results change in the DF approach. 




Fig. 2. Crossover temperature T* (mK) as a function of P (bar) of an 
electron bubble in liquid 4 He in the capillary model (dashed line) and in the 
DF approach (solid line). 



For a given temperature above T* , the solution to Eq. (jSJ) with Jt 
instead of Jq yields the pressure at which critical electron bubbles are created 
with sizeable probability. For the mentioned experimental parameters, it 
turns out that in the T ~ 1 K regime this happens when AU ~ 11.5 K. 
Inspection of Fig. ^ readily shows that this pressure is very close to P u . 
Consequently, P u is the key quantity for electron bubble explosionsP 

We would like to recall why the capillary model is expected to be fairly 
realistic for electron bubble explosions whereas it is completely unrealistic 



Cavitation of Electron Bubbles in Liquid Helium 



to address homogeneous cavitation in liquid helium. Dropping the elec- 
tron contributions to Eq. (^Q), it is easy to see that for negative pressures 
AU = 167t<t 3 /3|P| 2 at Rmax = 2a/\P\. This means that the barrier height is 
not zero at the spinodal pressure, which is very unphysical. This drawback 
renders useless the capillary model for homogeneous cavitation in pure liquid 
helium at low T because it has been experimentally established that cavita- 
tion occurs near the spinodal region EEE21 On the other hand, near R m i n the 
model yields empty bubbles of very small radius, which is also unphysical. 
As a consequence, any dynamical approach based on the use of these bubble 
configurations is rather dubious! 3 -^ The model can only work well for large 
bubbles, as for example near the saturation curve. In these situations, it has 
been shown that it yields realistic homogeneous cavitation pressures) 1 * 4 * 31 * 
and the use of Eq. © to obtain T* yields values in good agreement with 
the DF ones LLil2ll The situation is quite similar to that found in the descrip- 
tion of supersaturated 3 He- 4 He liquid mixtures, see Refs. ED and [331 for a 
thorough discussion. 



3. DENSITY-FUNCTIONAL APPROACH TO ELECTRON 

BUBBLE EXPLOSIONS 

Density functional theory has been applied in the past to cavitation 
in classical liquids (see Ref. and Refs. therein). Since the pioneering 
work of Xiong and Marisp'it has proven to be the most successful approach 
in addressing cavitation in liquid helium so far.^ It incorporates in a self- 
consistent way the equation of state of bulk liquid and surface tension of the 
liquid-gas interface as a function of temperature, which are key ingredients 
of any cavitation model. It does not impose a priori the density profile of 
the critical cavity, allowing for a flexible description of the process from the 
saturation line down to the spinodal line. Moreover, within DF theory one 
avoids to split the system into a bulk and a surface region, and the use of 
macroscopic concepts such as surface tension and pressure at a nanoscopic 
scale. However, it is a continuous, not an atomic description of the sys- 
tem. In spite of this, it has been found to well describe situations in which 
the atomic scale is relevant, such as quantized vortices, or the presence of 
strongly attractive atomic or molecular impurities (see for instance Ref. EH] 
and Refs. therein). 

In the frame of DF theory, the properties of an electron bubble ap- 
proaching the surface of liquid 4 He have been studied by Ancilotto and 
ToigoP^ They have used the so-called Orsay-Paris zero temperature finite- 
range DFp^and the pseudopotential proposed in Ref. [38 as e-He interaction. 
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The method chosen by Classen et alP is a simplification of that of Ref. EH1 
in two respects. Firstly, they have used a zero-range DF to describe 4 He, 
which seems justified in view of the slowly varying helium densities even 
in the presence of excess electrons in the liquid (the situation is completely 
different for strongly attractive atomic or molecular impurities, see for ex- 
ample Ref. I39j) . However, they have included some thermal effects in the 
DF, whereas the approach of Ancilotto and Toigo is at zero temperature.^ 
Secondly, the pseudopotential has been replaced by a contact e-He interac- 
tion whose intensity has been adjusted so as to reproduce the total energy 
of an excess electron in bulk helium,^ which is about 1 eV for a particle 
density of 0.0218 A~ 3 (saturation density). 

Our starting point is a finite temperature zero-range DF that reproduces 
thermal properties of liquid 4 He such as the experimental isotherms and the 
4 He liquid-gas coexistence line up to T = 4.5 K, and the T dependence of 
the surface tension of the liquid free surface.^ This DF has been successfully 
used to address homogeneous cavitation in liquid helium from T ~ K up to 
temperatures close to the critical one) 4 ^ In the quantum cavitation regime, it 
has also yielded results in good agreement with experiment and with other 
theoretical approaches.^ We have taken the Hartree-type e-He effective po- 
tential derived by Cheng et alP^l (see also Ref. I18|) as e-He interaction. This 
allows us to write the free energy of the system as a functional of the 4 He 
particle density p, the excess electron wave function \&, and T: 

F[o,-$,T] = I dff(p,T) + ^— [ cF|W(f)| 2 + / dr \^(r)\ 2 V(p) , (9) 
J 2 m e J J 

where f(p,T) is the 4 He free energy density per unit volume written as 

Hp, t) = f vo i(p, t) + + £(v P ) 2 . (io) 

p 

In this expression, f vo i(p,T) consists of the well-known free energy density 
of a Bose gas, plus phenomenological density dependent terms that take 
into account the effective interaction of helium atoms in the bulk liquid.^ 
The parameters of these terms and those of the density gradient terms in 
Eq. (fTU)) have been adjusted so as to reproduce physical quantities such 
as the equation of state of the bulk liquid and the surface tension of the 
liquid free surface. We have slightly modified the original value 2 -^ of the 
parameter £ in Eq. (|10|) to exactly reproduce the surface tension of liquid 
4 He (Ref. E0, taking £ = 2330 K A 5 . The /3-term is a kinetic energy term; 
at T = the system is described as a Bose condensate and for this reason 
the kinetic energy arises only from the inhomogeneity of the density.^ For 
inhomogeneous systems, this term is essential to have densities well behaved 
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everywhere, and in the case of 4 He droplets its inclusion in the functional 
yields densities that smoothly -exponentially- go from the bulk liquid down 
to zero. We have takerP/3= (ft 2 /2m 4 )/4. 

In bulk helium, knowledge of f(p, T) enables solution of the phase equi- 
librium equations and to determine the spinodal line. It also yields an 
equation of state in the negative pressure regime, inaccessible to the ex- 
perimental determination, through the thermodynamic relationship P = 
— fvoi{p,T) + fip, where \i is the 4 He chemical potential. 

The e-He interaction V(p) is written as a function of the local helium 
density^ 

V(p) = —^ + pa a -27rae 2 — p^ 3 , (11) 

2m e m e V o J 

where a = 0.208 A 3 is the static polar izability of a 4 He atom, and ko is 
determined from the helium local Wigner-Seitz radius r s = (3/47T/9) 1 / 3 by 
solving the transcendent equation 

tan[k (r s - a c )\ = k r s , (12) 

with a c and a a being the scattering lengths arising from a hard-core and 
from a polarization potential. We have takerPZlaQ, = —0.06 A, a c = 0.68 A. 

The application of DF theory to the cavitation problem proceeds as 
follows. For given P and T values one first determines the metastable and 
unstable cavities that would correspond to the local minimum and maximum 
configurations in the capillary model (actually, in the multidimensional space 
spanned by the more flexible DF configurations, the latter is no longer a 
local maximum but a saddle point). This is achieved by solving the Euler- 
Lagrange equations which result from the variation of the constrained grand 
potential density &(p, T) = u>(p, \P, T) — s\^\ 2 , where the grand potential 
density u(p,^, T) is defined from Eq. @ as 

u(p,V,T) = f(p,T) + ^—\W\ 2 + \y\ 2 V(p)-pp. (13) 



It yields 



Sf IT. 2^ 



2 m, 



A^ + V{p)y = e^, (15) 



where e is the lowest eigenvalue of the Schrodinger equation obeyed by the 
electron. These equations are solved assuming spherical symmetry, imposing 
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for p the physical conditions that p'(0) = and p{r — > oo) = pj,, where p b 
is the density of the metastable bulk liquid, and that the electron is in the 
Is state. Fixing p b and T amounts to fix P and T, as the pressure can be 
obtained from the bulk equation of state P = —f v d{Pb>T) + pp b , as well as 
p. Thus, p = df vo i(p,T)/dp\T is known in advance, whereas e is not and 
has to be determined from Eq. (|15jl. 
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Fig. 3. 4 He density profiles in A~ 3 (solid lines, right scale) and excess 
electron squared wave functions |^| 2 in A~ 3 (dashed lines, left scale), as a 
function of radial distance r (A). The top panel shows the stable bubble 
at T = K, P = bar. The other three panels show the near-to-unstable 
electron bubble for several (P, T) values. In the upper panel, the vertical thin 
dashed line indicates the radius of the capillary model bubble {R m in = 18.9 
A) , and the vertical thin solid line indicates the radius at which the helium 
density equals p b /2 \R{p b /2) = 18.0 A, with p b = 0.0218 A" 3 ]. 
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We have solved Eqs. (|14|) and (|15|) using a multidimensional Newton- 
Raphson method^ after having discretized them using n-point formulas for 
the r derivatives. We have used n = 13 formulas, but comparable results 
have been obtained using n = 7 and 9 formulas.^' A fine mesh of step 
Ar = 0.1 A has been employed, and the equations have been integrated up 
to i?oo = 150 A to make sure that the asymptotic bulk liquid has been 
reached. The multidimensional Newton-Raphson method has been applied 
until the local chemical potential -left hand side of Eq. (|14j) - does not differ 
substantially from p. We have checked that, for every r value, both coincide 
up to at least the sixth decimal figure. This is a crucial test on the accuracy 
of our method. We have thus achieved a fully variational solution of the 
Euler-Lagrange problem embodied in Eqs. (|14j) and (|15|) . valid from r = 
up to Roo. 




i i i i i i i i i i i 

1 2 3 4 5 

T(K) 

Fig. 4. Instability pressure P u as a function of T for 4 He. Circles and 
triangles are experimental data from Refs. |S]and|ni respectively. The results 
of the capillary model are represented by a dashed line, and the DF results 
by a solid line. The dash-dotted line represents the experimental saturation 
vapor pressure P sv line. The crosses at T = 2 and 4 K indicate the critical 
pressures P cr . 

We represent in Fig. several 4 He density profiles and excess electron 
squared wave functions l^] 2 . The top panel shows the stable bubble at T = 
K, P = bar. The other three panels display the near-to-unstable electron 
bubble for several (P, T) values. For a given T, they have been obtained 
decreasing pb, i.e. P, until Eqs. (jl"!]) and (fl3|) have no solution. The smaller 
P value defines P u . The 4 He instability pressure P u is shown in Fig. |l]as 
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a function of T. This figure shows that the lowest pressure the system may 
reach before becoming macroscopically unstable is P u = —2.07 bar, which is 
the value corresponding to T = K. 
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Fig. 5. 4 He density profiles in A 
tron squared wave functions 

of radial distance r (A) for the metastable and saddle electron bubble at 
different (P, T) values. Dashed lines correspond to the metastable configu- 
ration, and solid lines to the saddle configuration. From top to bottom, the 
^^max values are 65.3, 80.2, and 88.1 K, respectively. 



As we have indicated, in the metastability region the Euler-Lagrange 
equations have two different solutions for given pb and T -i.e., P and T- 
values, one corresponding to the metastable configuration and another cor- 
responding to the saddle configuration. A similar situation is found in the 
case of cavitation in the presence of vortices.™ Actually, it is the search of 
the saddle configuration that constitutes a numerical challenge. The reason 
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is that, due to the strongly repulsive e-He interaction, fairly small changes in 
the electron configuration induce changes in the helium configuration hard 
to handle numerically, so that in the course of the numerical procedure, the 
system has a strong tendency to jump from the saddle to the metastable 
solution. We have represented in Fig. [3 several 4 He density profiles and 
excess electron squared wave functions corresponding to the metastable 
and saddle electron bubbles for T = 0, 2 and 4 K, and a value of P close to 

The cavitation barrier height AQ max is obtained by subtracting the 
grand potential of the saddle bubble co(p s , *& S ,T) from that of the metastable 
bubble Lu(p m ,^ m ,T): 

AVt max = J df[oj(p s , # s , T) - u(p m , * m , T)] . (16) 

Since both configurations go asymptotically to the same pb value, Eq. 
(|16|) gives A£l max as a function of P and T via the equation of state of bulk 
liquid helium. AQ max is shown in Fig. 0as a function of P for T = 2 and 
4 K, and in Fig. □ for T = K. As indicated, An 

max becomes negligible 

when the system approaches the unstability pressure. This constitutes a 
suplementary test on the correctness of the near-to-unstable configurations 
we have found by decreasing pb at fixed T. Otherwise, AVL max would not be 
negligible for this configuration. 




P(bar) 

Fig. 6. DF energy barrier height AQ max (K) for He as a function of P 
(bar) for T = 2 K (left panel) and T = 4K (right panel). 
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Once AQ maxi-P^T) has been determined, it can be used to obtain the 
critical pressure P cr at which critical bubbles nucleate at an appreciable rate 
by solving an equation similar to Eq. 

l = t exp V exp J 0T e- A ^ p ^ T (17) 

taking Jqt = V T n e . This yields P cr = —1.55 bar at T = 2 K, and P cr = 0.211 
bar at T = 4 K. These pressures are slightly above the corresponding P u 
values which are, respectively, —1.62 and 0.151 bar. 

We thus see from Fig. HI that our results are fully compatible with the 
available experimental data -the calculated P u should be a lower bound to 
the cavitation pressure- if the experimental results of Ref. [Blare ruled out, 
and those carried out more recently by the same groupP have error bars 
similar to those found in the case of pure liquid 4 HeJ^ as well as in the case 
of electron bubble cavitation in 

3 He P 

The capillary model yields instability pressures lower than those ob- 
tained within DF theory, especially at high temperatures. The discrepancy 
arises even if one takes into account, as we did, the T-dependence of the 
surface tension. The origin of the discrepancy is the e-He strongly repulsive 
interation acting on the helium density tail that penetrates inside the cavity 
as T increases: the electron has to 'push' not only the bulk surface, which 
for DF configurations is at the radius where the helium density equals Pf,/2, 
but also the part of the helium density that is spread inside the bubble, 
especially at high temperatures. This 'pushing' produces a steeper bubble 
density profile, and not only an actual displacement of its surface. For this 
reason, the radius of the bubbles are smaller in the DF than in the simple 
capillary approach (see Figs. E]and|HI). Part of the difference is removed if 
one allows for penetration of the electron wave function into the liquid using 
a finite height barrier for the potential that confines the electron within the 
bubble. This would diminish the bubble radius, as it lowers the zero-point 
energy of the electron. However, the final effect is that it increases \P U \ and 
the agreement with the DF result worsens. ^ 

Our results for P u are slightly below those of Ref. OH (compare for in- 
stance our value of -2.07 bar at T = K with their value of -1.92 bar). 
This difference is essentially due to the different value of the surface tensions 
used to build the DF used in our work and in theirs. Indeed, the pressure 
we obtain if the DF is adjusted to reproduce the surface tension of Ref. |2U1 
is P u = —1.93 bar, in excellent agreement with the value found in Ref. 03 
Moreover, the ratio 2.07/1.93 compares very well with that obtained from 
Eq. © if ones takes o = 0.272 K A -2 in one case, and a = 0.257 K A~ 2 
in the other. We recall that in Ref. ED the kinetic /3-term in Eq. IflOj) has 
been neglected, which makes the helium density to be zero at the origin, 
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and that these authors have made the helium density strictly zero inside a 
sphere around the excess electron, whereas in our case the helium density is 
defined everywhere. These differences do not seem to play any role. 

We have also obtained AQ max at T = K using the contact e-He 
interaction used in Ref. 03 The results, indicated by crosses in Fig. 
indicate that the contact interaction sensibly yields the same barrier heights 
as the Hartree-type interaction. 

We have employed the DF theory to obtain the crossover temperature 
T*W^ To this end, one has to obtain the frequency of the small amplitude 
oscillations around the saddle configuration in the inverted barrier potential 
well. From the capillary model results we expect that T* is very small, so 
that thermal effects on p s can be neglected. In view of the fairly large size of 
the saddle bubble, the oscillation frequency can be obtained as follows E2I48I 
After determining the saddle configuration p s (r), we define a continuous set 
of densities by a rigid translation of /O s (r): 



The variable 5 roughly represents the displacement of the surface of the 
saddle configuration with respect to its stationary value. Varying 5, all 
physically relevant configurations are generated. Eq. (|18|) implies that the 
surface diffuseness of the bubble is kept frozen during the displacement. 
The barrier is then obtained as a function of 5: 



Within this model, 5 is the only collective variable describing the bubble 
oscillation, and all the time-dependence is in 5(t). The kinetic energy asso- 
ciated with the oscillation is 



where u(f, t) is the velocity field which can be formally obtained from the 
continuity equation 




(18) 




(19) 




(20) 



-£ + V(pu) = 0. 



(21) 



It yields: 




(22) 



By construction, 



p(r,t) =p s (r-5(t)). 



(23) 
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Thus 
and 



u(r, t) 



r 2 p s {r) 

Denning the mass parameter M(5) as 



p(r,t) = -p' s (r)S 



r 2 p s (r) 



sps(s)ds 



we get 



M(5) 



K 2 



dr 



r 2 p 5 (r) 



\m{5)5 2 



r 2 p s (r) 



2 / dssps(s) 



Proceeding as in the capillary case one obtains 



h 

2tt' 



d 2 AQ 



d5 2 



M{5) 



6=0 



(24) 
(25) 

(26) 

(27) 

(28) 



which is the generalization of Eq. (JSJ) to the case of diffuse density profiles. 
Using this equation we have obtained the T*{P) curve shown in Fig. |2 
and proceeding as in the capillary model, we have determined a crossover 
temperature of 6.0 mK. 



4. ELECTRON BUBBLE EXPLOSIONS IN LIQUID 3 He 

We have also studied the explosion of electron bubbles in the case of 
liquid 3 He. The capillary model of Sec. 2 can be straightfowardly applied to 
this isotope using the appropriate values of the surface tension a = 0.113 K 
A~ 2 (Ref. 1491 where one may also find the values of <r(T) we have used to 
obtain the capillary model results we show in Fig. EJ) and of the dielectric 
constant e = 1.0428 (Ref. 50) As in Sec. 2, we have neglected this term in 
the calculations because of its smallness. 

At zero pressure and temperature, the capillary model yields an electron 
bubble of radius R m i n = 23.5 A, larger than for 4 He because of the smaller 
3 He surface tension. For the same reason, the 3 He instability pressure P u = 
—0.71 bar is smaller in absolute value [see Eq. (j2J)]. This pressure is atteined 
for an electron bubble of radius R u = 35.2 A. 

The application of DF theory to describe electron bubble explosions in 
3 He proceeds as indicated in Sec. 3. We have used the functional proposed 
in Ref. 1511 which we have employed in the paslP^to describe homogeneous 
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Fig. 7. 3 He instability pressure P u as a function of T for 3 He. Dots are 
experimental data from Ref. 1101 The dash-dotted line represents the ex- 
perimental saturation vapor pressure P sv line. The results of the capillary 
model are represented by a dashed line, and the DF results by a solid line. 



cavitation in liquid 3 He, and the e-He interaction given in Eq. pij) with the 
parameters corresponding to 3 He, namely a = 0.206 A 3 , and same values 
for a a and a c . 

The 3 He instability pressure P u is shown in Fig. [7| as a function of 
T. This figure shows that the lowest pressure the system may reach before 
becoming macroscopically unstable is P u = —0.69 bar, which is again the 
value corresponding to T = K. For the reason indicated before, we have 
not considered necessary to calculate the critical pressures in the case of 3 He. 
It can be seen that our values of P u are somewhat above the experimental 
values of P cr obtained in Ref. 1101 especially at T ~ 1 K. As in the 4 He case, 
the capillary model fails as soon as thermal effects start being sizeable. 

We represent in Fig. El several 3 He density profiles and excess electron 
squared wave functions l^] 2 . The top panel shows the stable bubble at T = 
K, P = bar. The other three panels display the near-to-unstable electron 
bubble for several (P, T) values. Comparing with Fig. |31 it can be seen that 
3 He electron bubbles are more diffuse than 4 He electron bubbles. 

We have also obtained T* for 3 He in the capillary model. In spite of 
theoretical predictions that point to a crossover temperature of the order of 
100 mK in pure liquid 3 Hej 41 * 4 ^ recent experiments have not found such a 
crossing^, indicating that super fluid coherence might play a role in quantum 
cavitation. Yet, using the expressions given in Sec. 2 we have obtained 
T*(P) for 3 He and show it in Fig. |UJ It can be seen that the maximum of 
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T*(P) is roughly half the value we have obtained for 4 He in the capillary 
model. Proceeding as in the case of 4 He and assuming the same experimental 
parameteres, we have found that the value of T* that would correspond to 
this situation is about 2.3 mK. This estimate is close to the normal-to- 
superfluid transition temperature in 3 He, below which our DF method does 
not apply as it assumes that 3 He is in the normal phase. 
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Fig. 8. 3 He density profiles in A -3 (solid lines, right scale) and excess 
electron squared wave functions |^| 2 in A -3 (dashed lines, left scale), as a 
function of radial distance r (A). The top panel shows the stable bubble 
at T = K, P = bar. The other three panels show the near-to-unstable 
electron bubble for several (P, T) values. In the upper panel, the vertical thin 
dashed line indicates the radius of the capillary model bubble (R m in = 23.5 
A) , and the vertical thin solid line indicates the radius at which the helium 
density equals p b /2 [R(p b /2) = 22.5 A, with p b = 0.0163 A" 3 ]. 
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Fig. 9. Crossover temperature T* (mK) as a function of P (bar) of an 
electron bubble in liquid 3 He in the capillary model. 



5. SUMMARY 

We have thoroughly addressed electron bubble explosions in liquid he- 
lium. Our approach is based on the application of finite temperature den- 
sity functionals succesfully used to describe cavitation in liquid 4 He and 3 He. 
This approach is fully selfconsistent and unbiased by numerical artifacts, and 
to our knowledge, it is the only one applied in a wide range of temperatures. 

We have compared our results with experiments and have found that 
our calculations are either in agreement with the experimental datap^ or 
compatible with them if they have error bars similar to these attributed to 
other cavitation processes in liquid helium. 

We have used a realistic electron-helium interaction and have tested 
another approach based on the use of a simpler interaction. We have found 
that in spite of the fairly large electron bubbles involved in the process, 
at high temperatures the capillary model fails to yield quantitative results, 
overestimating the critical pressures. We attribute this to the 'rigidity' of the 
bubble configurations which is inherent to the capillary approach. Whereas 
these are serious drawbacks for nanoscopic bubbles, they are expected not 
to have a sizeable influence for microscopic multielectron bubbles 

We have also used the density functional results in conjunction with 
a functional integral method to obtain the thermal to quantum crossover 
temperature. This approach has led in the past to a correct description 
of the same process in pure liquid 4 He. In the present case, the crossover 
temperature turns out to be very small, about 6 mK. 
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Finally, we want to stress the suitability of the DF approach to quan- 
titatively address electron bubbles in liquid He. This might encourage one 
to investigate other problems like the infrared spectrum of the electron bub- 
ble in liquid helium, and the effect of quantized vortices pinned to excess 
electrons on the critical cavitation pressure. It has been arguecP^ that the 
rising of P u below T = 1 K could be attributed to the presence of quantized 
vortices. However, only simple models have been used to study their effect, 
and the agreement with experiment is only qualitative. In the case of the 
infrared spectrum of electron bubbles, the DF approach might shed light on 
the long-standing problem of how to understand the experimental results of 
Grimes and AdamJ^ on the Is — lp and Is — 2p electron transition ener- 
gies without using unjustifiable pressure dependences of the helium surface 
tension within the capillary model. 



ACKNOWLEDGMENTS 

We would like to thank F. Ancilotto, F. Caupin, M. Guilleumas, D. M. 
Jezek, H. Maris and F. Toigo for useful discussions. This work has been per- 
formed under grants BFM2002-01868 from DGI (Spain) and 2000SGR00024 
from Generalitat de Catalunya. 



REFERENCES 

1. J. A. Nissen, E. Bodegom, L. C. Brodie, and J. Semura, Phys. Rev. B 40, 617 
(1989). 

2. Q. Xiong and H. J. Maris, J. Low Temp. Phys. 82, 105 (1991). 

3. H. Lambare, P. Roche, S. Balibar, H. J. Maris, O. A. Andreeva, C. Guthmann, 
K. O. Keshishev, and E. Rolley, Eur. Phys. J. B 2, 381 (1998). 

4. S. Balibar, J. Low Temp. Phys. 129, 363 (2002). 

5. T. Satoh, M. Morishita, M. Ogata and S. Katoh, Phys. Rev. Lett. 69, 335 
(1992). 

6. S. Balibar, C. Guthmann, H. Lambare, P. Roche, E. Rolley, and H. J. Maris, 
J. Low Temp. Phys. 101, 271 (1995). 

7. X. Chavanne, S. Balibar, and F. Caupin, J. Low Temp. Phys. 126, 615 (2002). 

8. J. Classen, C.-K. Su, and H. J. Maris, Phys. Rev. Lett. 77, 2006 (1996). 

9. J. Classen, C.-K. Su, M. Mohazzab, and H. J. Maris, Phys. Rev. B 57, 3000 
(1998). 

10. C.-K. Su, C. E. Cramer, and H. J. Maris, J. Low Temp. Phys. 113, 479 (1998). 

11. F. Dalfovo, Phys. Rev. B 46, 5482 (1992). 

12. H. J. Maris, J. Low Temp. Phys. 94, 125 (1994). 

13. M. Barranco, M. Guilleumas, D. M. Jezek, R. J. Lombard, J. Navarro, and M. 
Pi, J. Low Temp. Phys. 117, 81 (1999). 



Cavitation of Electron Bubbles in Liquid Helium 



14. M. Barranco, M. Guilleumas, M. Pi, and D. M. Jezek, in Microscopic Ap- 
proaches to Quantum Liquids in Confined Geometries, E. Krotscheck and J. 
Navarro Eds. p. 319 (World Scientific, Singapore, 2002). 

15. M. Pi, M. Barranco, and R. Mayol, contribution to the Symposium on Quantum 
Fluids and Solids, Trcnto (Italy), July 5-9, 2004. J. of Low Temp. Phys., in 
print. Due to an unfortunate mistake in the electron-helium interaction energy, 
the numerical results presented in this reference are in error. 

16. R. A. Ferrel, Phys. Rev. 108, 167 (1957). 

17. D. Konstantinov and H. J. Maris, J. Low Temp. Phys. 121, 609 (2000). 

18. M. Rosenblit and J. Jortner, Phys. Rev. B 52, 17461 (1995). 

19. P. Roche, G. Deville, N. J. Appleyard, and F. I. B. Williams, J. Low Temp. 
Phys. 106, 565 (1997). 

20. M. lino, M. Suzuki, and A. J. Ikushima, J. Low Temp. Phys. 61, 155 (1985). 

21. A. Guirao, M. Centelles, M. Barranco, M. Pi, A. Polls, and X. Vinas, J. Phys. 
Condens. Matter A, 667 (1992). 

22. H. M. Guo, D. O. Edwards, R. E. Sarwinski, and J. T. Tough, Phys. Rev. Lett. 
27, 1259 (1971). 

23. Q. Xiong and H. J. Maris, J. Low Temp. Phys. 77, 347 (1989). 

24. J. Boronat, J. Casulleras, and J. Navarro, Phys. Rev. B 50, 3427 (1994). 

25. D. M. Jezek, M. Guilleumas, M. Pi, M. Barranco, and J. Navarro, Phys. Rev. 
5 48, 16 582 (1993). 

26. S. Coleman, Phys. Rev. D 15, 2929 (1977); C. G. Callan and S. Coleman, ibid. 
16, 1762 (1977). 

27. E. M. Chudnovsky, Phys. Rev. A 46, 8011 (1992). 

28. C.-K. Su and H. J. Maris, J. Low Temp. Phys. 110, 485 (1998). 

29. F. Caupin and S. Balibar, Phys. Rev. B 64, 064507 (2001). 

30. T. Nakamura, Y. Kanno, and S. Takagi, Phys. Rev. B 51, 8446 (1995). 

31. D. N. Sinha, J. S. Scmura, and L. C. Brodie, Phys. Rev. A 26, 1048 (1982). 

32. M. Guilleumas, M. Barranco, D. M. Jezek, R. J. Lombard, and M. Pi, Czech. 
J. Phys. 46, Suppl. SI, 389 (1996). 

33. E. Tanaka, K. Hatakeyama, S. Noma, S. N. Burmistrov, and T. Satoh, J. Low 
Temp. Phys. 127, 81 (2002). 

34. D. W. Oxtoby, J. Phys. Condens. Matter 4, 7627 (1992). 

35. F. Ancilotto, M. Barranco, and M. Pi, Phys. Rev. Lett. 91, 105302 (2003). 

36. F. Ancilotto and F. Toigo, Phys. Rev. B 50, 12 820 (1994). 

37. J. Dupont-Roc, M. Himbert, N. Pavloff, and J. Treiner, J. Low Temp. Phys. 
81, 31 (1990). 

38. N. R. Kestner, J. Jortner, M. H. Cohen, and S. A. Rice, Phys. Rev. 140, A56 
(1965). 

39. F. Dalfovo, Z. Phys. D 29, 61 (1994). 

40. W. T. Sommer, Phys. Rev. Lett. 12, 271 (1964); M. A. Woolf and G. W. 
Rayfield, ibid. 15, 235 (1965). 

41. H. J. Maris, J. Low Temp. Phys. 98, 403 (1995). 

42. E. Cheng, M. W. Cole, and M. H. Cohen, Phys. Rev. B 50, 1136 (1994); Erra- 
tum ibid. 50, 16134 (1994). 

43. S. Stringari and J. Treiner, J. Chem. Phys. 87, 5021 (1987); Phys. Rev. B 36, 
8369 (1987). 



M. Pi, M. Barranco, R. Mayol, and V. Grau 



44. W. H. Press, S. A. Teulosky, W. T. Vctterling, and B. P. Flannery, Numerical 
Recipes in Fortran 77: The Art of Scientific Computing (Cambridge University 
Press, Cambridge, 1999). 

45. W.G. Bickley, Math. Gaz. 25, 19 (1941); M. Abramowitz and I. A. Stegun, 
Handbook of Mathematical Functions (Dover, New York, 1965). 

46. D. M. Jezek, M. Guilleumas, M. Pi, and M. Barranco, Phys. Rev. B 55, 11092 
(1997). 

47. M. Guilleumas, M. Barranco, D. M. Jezek, R. J. Lombard, and M. Pi, Phys. 
Rev. B 55, 11092 (1997). 

48. M. Guilleumas, Ph. D. Thesis, University of Barcelona (1995), unpublished. 

49. M. lino, M. Suzuki, A. J. Ikushima, and Y. Okuda, J. Low Temp. Phys. 59, 
291 (1985). 

50. M. Rosenblit and J. Jortner, Phys. Rev. Lett. 75, 4079 (1995). 

51. M. Barranco, M. Pi, A. Polls, and X. Vinas, J. Low Temp. Phys. 80, 77 (1990). 

52. F. Caupin, S. Balibar, and H. J. Maris, J. Low Temp. Phys. 126, 91 (2002). 

53. J. Tcmpere, I. F. Silvcra, and J. T. Devreese, Phys. Rev. B 67, 035402 (2003). 

54. J. Classen, C.-K. Su, M. Mohazzab, and H. J. Maris, J. Low Temp. Phys. 110, 
431 (1998). 

55. C. C. Grimes and G. Adams, Phys. Rev. B 41, 6366 (1990). 



